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Abstract 

We consider the entanglement in a two dimensional XY model in an external magnetic field 
h. The model consists of a set of 7 localized spin-^ particles in a two dimensional triangular 
lattice coupled through nearest neighbor exchange interaction J. We examine the effect of single 
and double impurities in the system as well as the degree of anisotropy on the nearest neighbor 
entanglement and ergodicity of the system. We have found that the entanglement of the system 
at the different degrees of anisotropy mimics that of the one dimensional spin systems at the 
extremely small and large values of the parameter A = hj J . The entanglement of the Ising and 
partially anisotropic system show phase transition in the vicinity of A = 2 while the entanglement 
of the isotropic system suddenly vanishes there. Also we investigate the dynamic response of the 
system containing single and double impurities to an external exponential magnetic field at different 
degrees of anisotropy. We have demonstrated that the ergodicity of the system can be controlled 
by varying the strength and location of the impurities as well as the degree of anisotropy of the 
coupling. 
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I. INTRODUCTION 



Quantum entanglement is a corner stone in the structure of quantum theory with no 
classical analog Q . Entanglement is a nonlocal correlation between two (or more) quantum 
systems such that the description of their states has to be done with reference to each other 
even if they are spatially well separated. Particular fields where entanglement is considered 
as a crucial resource are quantum teleportation, cryptography and quantum computation 
, where it provides the physical basis for manipulating the linear superposition of the 
quantum states used to implement the different computational algorithms. On the other 
hand, many questions regarding the behavior of the complex quantum systems significantly 
rely on a deep understanding and a good quantification of the entanglement j4-9|. Partic- 
ularly, entanglement is considered as the physical resource responsible for the long range 
correlations taking place in many-body systems during quantum phase transitions. There 
has been great interest in studying the different sources of errors in quantum computing and 



their effect on quantum gate operations [lO|, lUl . Different approaches have been proposed 



for protecting quantum systems during the computational implementation of algorithms 
such as quantum error correction [l^ and decoherence-free subspace Q • Nevertheless, 
realizing a practical protection against the different types of induced decoherence is still 
a hard task. Therefore, studying the effect of naturally existing sources of errors such as 
impurities and lack of isotropy in coupling between the quantum systems implementing the 
quantum computing algorithms is a must. Furthermore, considerable efforts should be de- 
voted to utilizing such sources to tune the entanglement rather than eliminating them. The 



effect of impurities and anisotropy o: 
spin system has been investigated 



coupling between neighbor spins in a one dimensional 



151 . It was demonstrated that the entanglement can 



be tuned in a class of one-dimensional systems by varying the anisotropy of the coupling 
parameter as well as by introducing impurities into the spin system. For a physical quantity 
to be eligible for an equilibrium statistical mechanical description it has to be ergodic, which 
means that its time average coincides with its ensemble average. To test ergodicity for a 
physical quantity one has to compare the time evolution of its physical state to the corre- 
sponding equilibrium state. There has been an intensive efforts to investigate ergodicity in 
one-dimensional spin chains where it was demonstrated that the entanglement, magnetiza- 
tion, spin-spin correlation functions are non-ergodic in Ising and XY spin chains for finite 
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number of spins as well as at the thermodynamic limit js], [igIiSI. 

Studying quantum entanglement in two-dimensional systems face more obstacles in com- 
parison to the one dimensional case, particularly the rapid increase in the dimension of the 
Hilbert spaces which lead to much larger scale calculations relying mainly on the numerical 
methods. The existence of exact solutions has contributed enormously to the understanding 
of the entanglement for ID systems a previous work, the entanglement in a 

19-site two-dimensional transverse Ising model at zero temperature j21| was studied. The 
spin- 1/2 particles are coupled through an exchange interaction J and subject to an external 
time-independent magnetic field h. It was demonstrated that for such a class of systems 
the entanglement can be tuned by varying the parameter \ = h/J and also by introducing 
impurities into the system, which showed a quantum phase transition at a critical value of 
the parameter A in the vicinity of 2. Recently, we have investigated the time evolution of 
entanglement in a two dimensional triangular transverse Ising system with seven spins in 
an external magnetic field j22|. Different time dependent forms of the magnetic field were 
applied. The system have demonstrated different responses based on the type of apphed 
field, where for a smoothly changing magnetic field the system entanglement follows the 
profile of the field very closely. 

In this paper, we consider the entanglement in a two-dimensional XFtriangular spin 
system, where the nearest neighbor spins are coupled through an exchange interaction J 
and subject to an external magnetic field h. We consider the system at different degrees of 
anisotropy to test its effect on the system entanglement and dynamics. The number of spins 
in the system is 7 with a number of impurities existing. We consider two different cases of 
impurities, the first case is a single impurity existing either at the border of the system or at 
the center with the coupling strength between the impurity spin and its neighbors different 
from that between the rest of the spins. The second case is double impurities, existing both 
at the border or one at the border and one at the center. We consider the coupling between 
the two impurities as J', which is different from the coupling J" between each one of them 
and its neighbors, while the interaction among the other spins is J. 

We show that the entanglement profile of the system at different degrees of anisotropy 
has great resemblance to that of the one dimensional spin systems as the parameter A — )■ 
and oo. On the other hand, both the Ising and the partially anisotropic systems show 
phase transition behavior in the vicinity of A = 2 but the isotropic system show sharp 
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step variations in the same region before suddenly vanishing. Examining the effect of an 
external exponential magnetic field on the time evolution of the entanglement showed that 
the ergodicity of the system can be tuned by varying the strength and location of the 
impurities and the degree of anisotropy in the system. 

This paper is organized as follows. In the next section we present our model and quan- 
tification of entanglement. In sec. Ill we consider the case of a single impurity. In sec IV 
we study the system with a double impurity. We conclude in sec V. 

II. MODEL AND QUANTIFICATION OF ENTANGLEMENT 

We consider a set of 7 localized spin-i particles in a two dimensional triangular lattice 
coupled through exchange interaction J and subject to an external time-dependent magnetic 
field of strength h{t). All the particles are identical except one (or two) of them which are 
considered impurities. The Hamiltonian for such a system is given by 
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FIG. 1: The two dimensional triangular spin lattice in presence of an external transverse magnetic field. 

except the sites nearest to an impurity site. For a single impurity, the coupling between 
the impurity and its neighbors Ji^j — J' — (a -|- 1) J, where a measures the strength of 
the impurity. For double impurities Jj^ — J' — (cti -|- 1) J is the coupling between the two 
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impurities and Jij = J" = (02 + 1)J is the couphng between any one of the two impurities 
and its neighbors while the coupling is just J between the rest of the spins. 

For this model it is convenient to set J = 1. For a system of 7 spins, its Hilbert space is 
huge with 2'' dimensions, yet it is exactly diagonalizable using the standard computational 
techniques. Exactly solving Schrodinger equation of the Hamiltonian ([T]), yielding the system 
energy eigenvalues Ei and eigenf unctions ipi. The density matrix of the system is defined by 

p=\^o){ipo\, (2) 
where \ipo) is the ground state energy of the entire spin system. We confine our interest to 



the entanglement between two spins, at any sites i and j [23[ . All the information needed in 
this case, at any moment t, is contained in the reduced density matrix Pij{t) which can be 
obtained from the entire system density matrix by integrating out all the spins states except 
i and j. We adopt the entanglement of formation, as a well known measure of entanglement 
where Wootters has shown that, for a pair of binary qubits, the concurrence C, which 
goes from to 1, can be taken as a measure of entanglement. The concurrence between two 
sites i and j is defined as 

C(p) = max{0, ei - £2 - £3 - 64}, (3) 



where the e^'s are the eigenvalues of the Hermitian matrix R = a/ \^Py/p with p = {a^ ® 
(j^)p*((7^ ® cr^) and is the Pauli matrix of the spin in y direction. For a pair of qubits the 
entanglement can be written as, 

E{p) = e{C{p)), (4) 
where e is a function of the "concurrence" C 

.,0 . . f , (5) 



where h is the binary entropy function 

h{x) = —X log2 X — (1 — x)log2{l — x). (6) 

In this case, the entanglement of formation is given in terms of another entanglement 
measure, the concurrence C. The dynamics of entanglement is evaluated using the same 



techniques applied in our previous work [22|. Specifically, we apply the step-by-step time- 



evolution projection technique, which was proved to give the same exact result as the matrix 



transformation technique, where both techniques were introduced in [22|, but 20 times faster. 
In this technique we assume that our system is initially, at to, in the ground state at zero 
temperature |0) with energy, say, e in an external magnetic field with strength a. The 
magnetic field is turned to a new value b and the system Hamiltonian becomes H with 
eigenpairs Ei and The original state |0) can be expanded in the basis 

10) = Cil^i) + C2|^2) + - + CnIiPn), (7) 

where 

C^ = m(P). (8) 

When H is independent of time between t and to then we can write 

U{t, to) = e-^(*>*«)(*-*»)/'"1^,,„) = e-^>(*-*«)/'"1^,,„), (9) 

where U {t, to) is the time evolution operator. The ground state will evolve with time as 

N 

= 5^Q|V^,)e-*^»(*-*«). (10) 
1=1 

and the pure state density matrix becomes 

p{t) = \m{m\- (11) 

Simply any complicated function can be treated as a collection of step functions. When the 
state evolves to the next step just repeat the procedure to get the next step results. Of course 
the lack of smoothness in the magnetic field function imposes a challenging obstacle in the 
calculations but this can be overcome by choosing a proper small enough time step. Because 
the size of our 7-site system is still manageable, in our actual calculations, we included all 
the 2^ = 128 states in every step, without any truncation of the higher energy eigenstates. 
This ensures us no approximation in this step. But the method itself is aiming at larger size 
system, like 19 sites XY model. By then, due to the computation limit, cutting off higher 
energy eigenstates might be a necessary action. 
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III. SINGLE IMPURITY 



A. Static system with border impurity 

We define a dimensionless coupling parameter \ = h/J and we set J = 1 throughout 
this paper for convenience. We start by considering the effect of a single impurity located 
at the border site 1. The concurrence between the impurity site 1 and site 2, C(l, 2), versus 
the parameter A for the three different models, Ising (7 = 1), partially anisotropic (7 = 0.5) 
and isotropic XY (7 = 0) at different impurity strengths (a = —0.5, 0, 0.5, 1) is in fig. |2l 
Firstly, the impurity paramter a is set to zero. For the corresponding Ising model, the 
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FIG. 2: (Color online) The concurrence C(l,2) versus the parameter A with a single impurity at the 
border site 1 with different impurity coupling strengths a — —0.5, 0, 0.5, 1 for different degrees of anisotropy 
7 = 1, 0.5, as shown in the subfigures. The legend for all subfigures is as shown in subfigure (a). 

concurrence C(l,2), in fig. |2]^a), demonstrates the usual phase transition behavior where it 
starts at zero value and increases gradually as A increases reaching a maximum at A ~ 2 



then decays as A increases further. As the degree of anisotropy decreases the behavior of 
the entanglement changes, where it starts with a finite value at A = and then shows a step 
profile for the small values of A. For the partially anisotropic case, the step profile is smooth 
and the entanglement mimics the Ising case as A increases but with smaller magnitude. The 
entanglement of isotropic XY system shows a sharp step behavior then suddenly vanishes 
before reaching A = 2. Interestingly, the entanglement behavior of the two-dimensional spin 
system at the different degrees of anisotropy mimics the behavior of the one-dimensional 
spin system at the same degrees of anisotropy at the extreme values of the parameter A. 
The ground state of the one- dimensional Ising model is characterized by a quantum phase 
transition that takes place at the critical value h/J = 1 8] which corresponds to a 
maximum entanglement in the system. The order parameter is the magnetization (cx^) 
which is different from zero for J > h and zero otherwise. The ground state is paramagnetic 
when J//i — > where the spins get aligned in the magnetic field direction, the ^-direction. It 
is ferromagnetic when J//i — )■ oo where the spins are aligned in the x-direction. Both cases 
cause zero entanglement. Comparing the entanglement behavior in the two-dimensional 
Ising spin system with the one-dimensional system, one can see a great resemblance except 
that the critical value becomes h/J^2 in the two dimensional case as shown in fig. [2l On 
the other hand, for the partially anisotropic and isotropic XY systems, the entanglements 
of the two-dimensional and one-dimensional system agrees at the extreme values of A where 
it vanishes ioi h » J and reaches a finite value for h « J . The former case corresponds 
to an alignment of the spins in the 2;-direction, paramagnetic state, while the latter case 
corresponds to alignment in the x and y-directions which a ferromagnetic state. 

The effect of a weak impurity (J' < J), a = —0.5, is shown in fig. [2](b) where the 
entanglement behavior is the same as before except that the entanglement magnitude is 
reduced compared with the pure case. On the other hand, considering the effect of a strong 
impurity (J' > J), where a = 0.5 and 1, as shown in fig. [2](c) and fig. [2](d) respectively, one 
can see that the entanglement profile for 7 = 1 and 0.5 have the same overall behavior as in 
the pure and weak impurity cases except that the entanglement magnitude becomes higher 
as the impurity gets stronger and the peaks shift toward higher A values. Nevertheless, the 
isotropic XY system behaves differently from the previous cases where it starts to increase 
first in a step profile before suddenly dropping to zero again, which will be explained latter. 
To study the entanglement between two sites, none of them is impurity, we consider C(2,4) 
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FIG. 3: (Color online) The concurrence C(2,4) versus the parameter A with a single impurity at the 
border site 1 with different impurity coupling strengths a — —0.5, 0, 0.5, 1 for different degrees of anisotropy 
7 = f , 0.5, as shown in the subfigures. The legend for all subfigures is as shown in subfigure (a). 

which is depicted in fig. |3l There are two main differences between the behavior of C(2,4) 
and C(l,2). Firstly, the magnitude of the entanglement envelope is higher for C(2,4) for 
7 = 0.5 and (but not 7 = 1) when a = 0, while C(2,4) is greater than C(l,2) for all 
7 values for the weak impurity case, a = 0. This is an interesting result as internal sites 
entanglement should be smaller in value than the edge sites. Secondly, the entanglement 
of the isotropic XY case increases in a multi-step profile for all values of a before suddenly 
dropping to zero. 
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B. Static system with center impurity 



To explore the effect of the impurity location we investigate the case of a single impurity 
spin located at site 4, instead of site 2, where we plot the concurrences C(l, 2) and C(l, 4) in 
figs, m and E] respectively. Interestingly, while changing the impurity location has almost no 
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FIG. 4: (Color online) The concurrence C(l,2) versus the parameter A with a single impurity at the 
central site 4 with different impurity coupling strengths a — —0.5, 0, 0.5, 1 for different degrees of anisotropy 
7 = 1, 0.5, as shown in the subfigures. The legend for all subfigures is as shown in subfigure (a). 



effect on the behavior of the entanglement C(l, 2) of the partially anisotropic and isotropic 
XY systems, it has a great impact on that of the Ising system where the peak value of the 
entanglement increases significantly in the weak impurity case and decreases as the impurity 
gets stronger as shown in fig. HI Now considering the entanglement between the central im- 
purity site 4 and the edge site 1, and comparing with the results in fig. [3] of the entanglement 
between the edge site 2 and central site 4, one can see that the entanglement C(l,4) profile 
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FIG. 5: (Color online) The concurrence C(l,4) versus the parameter A with a single impurity at the 
border site 4 with different impurity coupling strengths a — —0.5, 0, 0.5, 1 for different degrees of anisotropy 
7 = 1, 0.5, as shown in the subfigures. The legend for all subfigures is as shown in subfigure (a). 

for all degrees of anisotropy is very close to the C(2,4). Nevertheless the entanglement 
C(l, 4) magnitude is lower for weak impurity case and higher for the strong impurity which 
means that the central impurity made a significant change to the entanglement magnitude. 



C. Effect of system energy gap on entanglement 

To explain the distinct behavior of the entanglement corresponding to the different degrees 
of anisotropy 7 we depict the lowest few energy eigenvalues of the system at the different 7 
values for the two cases of border and central impurities in fig. [6] and fig. [7] respectively. As 
can be noticed in fig. EJ^a), the energies of the ground state and the first excited state of the 
Ising system coincide at the beginning at the small values of A until a specific value where 
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FIG. 6: (Color online) The energy spectrum versus the parameter A with a single impurity at the border 
site 1 with impurity coupling strength a = 1 for different degrees of anisotropy 7 = 1, 0.5, as shown in the 
subfigures. The legend for all subfigures is as shown in subfigure (a). 

they deviate from each other. This is corresponding to the transition from the degenerated 
ground state to non-degenerated one, from paramagnetic to ferromagnetic order, by breaking 
the Z2 symmetry, which explains the phase transition curve observed in the Ising case. The 
energy spectrum of the partially anisotropic XY system is a little bit different at the small 
values of A where the ground state and the first excited state coincide at the beginning 
but then deviates slightly from each other before recombining again and at last separate 
from each other completely, this behavior is repeated quite few times depending on the 
impurity strength, as illustrated in fig. Mjo). This energy spectrum behavior explains the 
roughness in the ascending part of the entanglement curves of the partially anisotropic 
XY system corresponding to subsequent transitions between the ground state and the first 
excited state taking place before reaching the maximum entanglement point. In fig. EJ^c), 
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FIG. 7: (Color online) The energy spectrum versus the parameter A with a single impurity at the central 
site 4 with impurity coupling strength a — I for different degrees of anisotropy 7 = 1,0.5,0 as shown in the 
subfigures. The legend for all subfigures is as shown in subfigure (a). 

the energy spectrum of the isotropic XY system is explored where clearly the deviations 
and recombination between the ground and first excited state energies become sharper and 
more frequent compared with the partially anisotropic system. This distinct behavior of the 
energy spectrum corresponding to 7 = is the reason for the sharp step behavior of the 
entanglement as was shown in figs. [2] and El Critical quantum behavior in a many body 
system happens either when an actual crossing takes place between the excited state and 
the ground state or a limiting avoided level-crossing between them exists, i.e. an energy gab 
between the two states that vanishes in the infinite system size limit at the critical point 
20l |. When a many body system crosses a critical point, significant changes in both its wave 
function and ground state energy takes place, which are manifested in the behavior of the 
entanglement function. The entanglement in one dimensional infinite spin systems, Ising 
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FIG. 8: (Color online) (a) The concurrence C14 versus A; (b) the first derivative of the concurrence C14 
with respect to A versus A; (c) the energy gap between the ground state and first excited state versus A; (d) 
the first derivative (in units of J) and second derivative (in units of J^) of the energy gap with respect A 
versus A for the pure Ising system (7 = 1 and a = 0). 



and XY ^ was shown to demonstrate scaling behavior in the vicinity of critical points |23 |. 
The change in the entanglement across the critical point was quantified by considering the 
derivative of the concurrence with respect to the parameter A. This derivative was explored 
versus A for different system sizes and although it didn't show divergence for finite system 
sizes, it showed clear anomalies which developed to a singularity at the thermodynamic limit. 
The ground state of the Heisenberg spin model is known to have a double degeneracy for 
an odd number of spins which is never achieved unless the thermodynamic limit is reached 
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FIG. 9: (Color online) (a) The concurrence C14 versus A; (b) the first derivative of the concurrence C14 
with respect to A versus A; (c) the energy gap between the ground state and first excited state versus A; (d) 
the first derivative (in units of J) and second derivative (in units of J^) of the energy gap with respect A 
versus A for the pure partially anisotropic system (7 — 0.5 and a — 0). Notice that the first derivative of 
energy gap is enlarged 10 times its actual scale for clearness. 



20| . Particularly, the Ising ID spin chain in an external transverse magnetic field has 



doubly degenerate ground state in a ferromagnetic phase that is gapped from the excitation 
spectrum by 2J(1 — h/J), which is removed at the critical point and the system becomes 
in a paramagnetic phase. Now let us first consider our two-dimensional finite size Ising spin 
system. The concurrence Cu and its first derivative are depicted versus A in figs. M^a) and 
(b) respectively. As one can see, the derivative of the concurrence shows strong tendency 
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FIG. 10: (Color online) (a) The concurrence C14 versus A; (b) the first derivative of the concurrence C14 
with respect to A versus A; (c) the energy gap between the ground state and first excited state versus A; (d) 
the first derivative (in units of J) and second derivative (in units of J^) of the energy gap with respect A 
versus A for the pure isotropic system (7 = and a = 0). Notice that the first derivative of energy gap is 
enlarged 20 times its actual scale for clearness. 

of being singular at Ac = 1.64. The characteristics of the energy gap between the ground 
state and the first excited state as a function of A are explored in fig. lU^c). The system 
shows strict double degeneracy, zero energy gap, only at A = i.e. at zero magnetic field, 
but once the magnetic field is on the degeneracy is lifted and an extremely small energy gap 
develops, which increase very slowly for small magnetic field values but increases abruptly 
at certain A value. It is important to emphasis here that at A = 0, regardless of which one 
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of the double ground states is selected for evaluating the entanglement, the same value is 
obtained. The critical point of a phase transition should be characterized by a singularity in 
the ground state energy, and an abrupt change in the energy gap of the system as a function 
of the system parameter as it crosses the critical point. To better understand the behavior 
of the energy gap across the prospective critical point and identify it, we plot the first and 
second derivatives of the energy gap as a function of A in fig. [8]^d). Interestingly, the first 
derivative dAE/dX which represents the rate of change of the energy gap as a function of A 
starts with a zero value at A = and then increase very slowly before it shows a great rate 
of change and finally reaches a saturation value. This behavior is best represented by the 
second derivative d'^AE/dX'^, which shows strong tendency of being singular at Ac = 1.8, 
which indicates the highest rate of change the energy gap as a function of A. The reason 
for the small discrepancy between the two values of the Ac extracted from the dC/dX plot 
and the one of d'^AE/dX^ is that the concurrence C14 is only between two sites and does not 
represent the whole system in contrary to the energy gap. One can conclude that the rate 
of change of the energy gap as a function of the system parameter, A in our case, should 
be maximum across the critical point. Turning to the case of the partially anisotropic spin 
system, 7 = 0.5, presented in fig. [9l one can notice from fig. [9t^a) that the concurrence 
shows few sharp changes, which is reflected in the energy gap plot as an equal number of 
minima as shown in fig. [9t^b). Nevertheless, again there is only one strict double degeneracy 
at A = while the other three energy gap minima are non-zero and in the order of 10^^. It 
is interesting to notice that the anomalies in both dC/dX and d'^AE/dX^ are much stronger 
and sharper compared with the Ising case as shown in figs.[9t^c) and (d). Finally the isotropic 
system which is depicted in fig. [TOl shows even sharper energy gap changes as a result of the 
sharp changes in the concurrence and the anomalies in the derivatives dC/dX and d'^AE/dX^ 
are even much stronger than the previous two cases. 

D. System dynamics with impurity 

Now we turn to the dynamics of the two dimensional spin system under the effect of a 
single impurity and different degrees of anisotropy. We investigate the dynamical reaction 
of the system to an applied time- dependent magnetic field with exponential form h{t) = 
b+{a- b)e-'" * for t > and h{t) = a for t < 0. 
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FIG. 11: (Color online) Dynamics of the concurrences C(l, 2), C(l, 4), C(2, 4) with a single impurity at the 
border site 1 with different impurity coupling strengths a — —0.5,0, 1,2 for the two dimensional Ising lattice 
(7 = 1) under the effect of an exponential magnetic field with parameters values a=l, b=3.5 and uj — 0.1. 
The straight lines represent the equilibrium concurrences corresponding to constant magnetic field h = 3.5. 
The legend for all subfigures is as shown in subfigure (a). 

We start by considering the Ising system, 7 = 1 with a single impurity at the border site 
1, which is explored in fig. [TT], where we set a = 1, 6 = 3.5 and u = 0.1. For the pure case, 
a = shown in fig. [TlTa). the results confirms the ergodic behavior of the system that was 
demonstrated in our previous work 22], where the asymptotic value of the entanglement 
coincide with the equilibrium state value at h{t) = b. As can be noticed from figs. [TlTb). 
[TlTc) and [TlT d) neither the weak nor strong impurities have effect on the ergodicity of the 
Ising system. Nevertheless, there is a clear effect on the asymptotic value of entanglements 
C(l,2) and C(l,4) but not on C(2,4) which relates two regular sites. The weak impurity, 
a = —0.5 reduces the asymptotic value of C(l,2) and C(l,4) while the strong impurities. 
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a = 1,2 raise it compared to the pure case. In fig. [12], we consider the same system but 
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FIG. 12: (Color online) Dynamics of the concurrences C(l, 2), C(l, 4), C(2, 4) with a single impurity at the 
border site 1 with different impurity coupling strengths a = —0.5,0, f,2 for the two dimensional Ising lattice 
(7 = 1) under the effect of an exponential magnetic field with parameters values a=f , b=f .5 and cj = 0.1. 
The straight lines represent the equilibrium concurrences corresponding to constant magnetic field h ~ f .5. 
The legend for all subfigures is as shown in subfigure (a). 



under the effect of a weaker exponential magnetic field with set of parameters a = 1, 6 = 1.5 
and = 0.1. As can be noticed, the entanglement of the Ising system is still showing 
an ergodic behavior at all impurity strengths. This means that the Ising system with a 
single border impurity is always ergodic under the effect of different impurity strengths 
and different magnetic field parameters. Interestingly, the different impurity strengths have 
different effects on the asymptotic value of the entanglements compared to the previous case 
under the effect of the new magnetic field. The weak impurity, as shown in fig. [T2](b). raises 
the asymptotic value of C(2,4) and splits those of C(l,2) and C(l,4) from each other. 
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On the other hand, the strong impurity effects are depicted in fig. [TW c) and fig. [TW d) 
which show that the asymptotic values of all concurrences are reduced significantly as the 
impurity strength increases, in contrary to the previous case. This emphasis the important 
role that the magnetic filed parameters play beside the impurity strength in controlling the 
entanglement behavior. It is of great interest to examine the effect of the impurity location. 
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FIG. 13: (Color online) Dynamics of the concurrences C(l, 2), C(l, 4), C(l, 5) with a single impurity at the 
central site 4 with different impurity coupling strengths a = — 0.5,0, 1,2 for the two dimensional Ising lattice 
(7 — 1) under the effect of an exponential magnetic field with parameters values a=l, b=:1.5 and w = 0.1. 
The straight lines represent the equilibrium concurrences corresponding to constant magnetic field h — 1.5. 
The legend for all subfigures is as shown in subfigure (a). 

which we investigate in fig. where a single impurity is located at the central site 4 instead 
of the border site 1 with exponential magnetic field parameters a = 1, 6 = 1.5 and u; = 0.1. 
Very interestingly, the entanglement behavior changes significantly as a result of changing 
the impurity location. Though the pure Ising system is still ergodic as shown in fig. [TST a). the 
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system with weak and strong impurity becomes non-ergodic which is illustrated in figs. [TST b). 
[T3](c) and [T3](d) respectively. Again the weak impurity raises the asymptotic values wheres 
the strong impurities reduces them significantly. The dynamics of the partially anisotropic 
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FIG. 14: (Color online) Dynamics of the concurrences C(l, 2), C(l, 4), C(2, 4) with a single impurity at the 
border site 1 with different impurity coupling strengths a — —0.5,0, 1,2 for the two dimensional partially 
anisotropic lattice (7 = 0.5) under the effect of an exponential magnetic field with parameters values a=l, 
b=3.5 and uj = 0.1. The straight lines represent the equilibrium concurrences corresponding to constant 
magnetic field h = 3.5. The legend for all subfigures is as shown in subfigure (b). 



XY system under the effect exponential magnetic field with parameters a = 1, b = 3.5 and 
uj = 0.1, is explored in fig. [141 It is remarkable to see that while for both the pure and 
weak impurity cases, a = and —0.5 , the system is nonergodic as shown in figs. [T^a) and 
fl^b). and it is ergodic in the strong impurity cases a = 1 and 2 as illustrated in figs. fWc) 
and [T^d). Clearly the asymptotic values of the entanglement are higher in the pure 
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FIG. 15: (Color online) Dynamics of the concurrences C(l, 2), C(l, 4), C(2, 4) with a single impurity at the 
border site 1 with different impurity coupling strengths a — —0.5,0, 1,2 for the two dimensional partially 
anisotropic lattice (7 = 0.5) under the effect of an exponential magnetic field with parameters values a=l, 
b=1.5 and cj = 0.1. The straight lines represent the equilibrium concurrences corresponding to constant 
magnetic field h — 1.5. The legend for all subfigures is as shown in subfigure (a). 

and weak impurity cases compared with the strong impurities. Changing the magnetic field 
parameter value b to 1.5, one can observe the great impact in fig. [151 where only the pure 
system becomes ergodic while the system with any impurity strength is nonergodic. This 
means that the magnetic field parameters control ergodicity as well. In fig. [161 we study 
the same system with a single impurity at the central site 4 with magnetic field parameters 
a = 1, = 1.5 and w = 0.1. As can be seen from the different subfigures, the system is 
ergodic only in the pure case and the asymptotic values are higher for the pure and weak 
impurity systems compared with the strong impurities. The complete isotropic XY system 
with a single border impurity at site 1 under the effect of an exponential magnetic field 
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FIG. 16: (Color online) Dynamics of the concurrences C(l, 2), C(l, 4), C(2, 4) with a single impurity at the 
central site 4 with different impurity coupling strengths a — —0.5,0,1,2 for the two dimensional partially 
anisotropic lattice (7 — 0.5) under the effect of an exponential magnetic field with parameters values 
a = 1,6 = 1.5 and u = 0.1. The straight lines represent the equilibrium concurrences corresponding to 
constant magnetic field h — 1.5. The legend for all subfigures is as shown in subfigure (a). 

with parameter values a = 1, 6 = 1.5 and u = 0.1, is investigated in fig. [T71 The trivial 
effect of the magnetic field, similar to the one dimensional case results |8|, is clear where 
the entanglement assumes a constant value for all pair of spins. This trivial effect is the 
result of the fact that for 7 = the exchange coupling terms in the Hamiltonian commute 
with the magnetic field term. Nevertheless one still can see an effect of the impurity on the 
ergodicity of the system where for a = and 1, the system is nonergodic wile for a = —0.5 
and 2 it is ergodic as shown in fig. [T71 In fact, testing a wide range of a values indicates 
that for the values approximately in the range —0.4 > a < 1.9 the system is nonergodic, 
otherwise it is ergodic i.e. for small absolute values of the impurity. Examining the same 
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FIG. 17: (Color online) Dynamics of the concurrences C(l, 2), C(l, 4), C(2, 4) with a single impurity at the 
border site 1 with different impurity coupling strengths a — —0.5, 0, 1, 2 for the two dimensional XY lattice 
(7 = 0) under the effect of an exponential magnetic field with parameters values a = 1, 6 = 1.5 and uj — 0.1. 
The straight (thicker) lines represent the equilibrium concurrences corresponding to constant magnetic field 
h — 1.5. The legend for all subfigures is as shown in subfigure (a). 

system under the effect of the same magnetic field but with a single central impurity, for 
wide range of a, demonstrates that the system becomes nonergodic at all values of a which 
is illustrated in fig. [TSJ 



24 



0.4 



0)0-3 

o 

c 

§ 0.2 
o 

o 
O 



0.1 



(a) 



20 



a = 



40 



60 



t (units of J"^) 



80 



100 



0.4 



CD 0-3 
o 

§ 0.2 
o 

o 

"0.1 



(b) 



20 



a = - 0.5 



40 



60 



t (units of J~^) 



-C(1,2) 
-C(1,4) 
C(1,5) 



80 



100 





0.5 


CD 




O 


0.4 




CD 
i_ 




i_ 


0.3 


o 




d 




o 
O 


0.2 




0.1 



(c) 



20 



a= 1 



40 



60 



t (units of J"') 



80 



100 



CD 



0.4 



0.3 



^ 0.2 
o 

o 

"o.i 



(d) 



20 



a = 2 



40 



60 



t (units of J~') 



80 



100 



FIG. 18: (Color online) Dynamics of the concurrences C(l, 2), C(l, 4), C(2, 4) with a single impurity at the 
central site 4 with different impurity coupling strengths a — — 0.5,0, 1,2 for the two dimensional XY lattice 
(7 = 0) under the effect of an exponential magnetic field with parameters values a = 1, 6 = 1.5 and uj — 0.1. 
The straight (thicker) lines represent the equilibrium concurrences corresponding to constant magnetic field 
h = 1.5. The legend for all subfigures is as shown in subfigure (b). 

IV. DOUBLE IMPURITIES 

A. Static system with impurities 



In this section we study the effect of double impurity, where we start with two located 
at the border sites 1 and 2. We set the coupling strength between the two impurities as 
J' = (1 + «i)J, between any one of the impurities and its regular nearest neighbors as 
J" = (1 + 02) J and between the rest of the nearest neighbor sites on the lattice as J. The 
effect of the impurities strength on the concurrence between different pairs of sites for the 
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FIG. 19: (Color online) The concurrence C(l, 2), C(l, 4), C(4, 5) versus the impurity coupling strengths ai 
and a2 with double impurities at sites 1 and 2 for the two dimensional Ising lattice (7 = 1) in an external 
magnetic field h=2. 

Ising lattice is shown in fig. [191 In fig. [TOT a) we consider the entanglement between the two 
impurity sites 1 and 2 under a constant external magnetic field h = 2. The concurrence 
C(l, 2) takes a large value when the impurity strengths ai, controlling the coupling between 
the impurity sites, is large and when 02, controlling coupling between impurities and their 
nearest neighbors, is weak. As ai decreases and ^2 increases, C(l, 2) decreases monotonically 
until it vanishes. As one can conclude, ai is more effective than a2 in controlling the 
entanglement in this case. On the other hand, the entanglement between the impurity 
site 1 and the regular central site 4 is illustrated in fig. [TW b) which behaves completely 
different from C(l,2). The concurrence C(l,4) is mainly controlled by the impurity strength 
0^2 where it starts with a very small value when the impurity is very weak and increases 
monotonically until it reaches a maximum value at 02 = 0, i.e. with no impurity, and 
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FIG. 20: (Color online) The concurrence C(l,2), C(l,4), C(2,5), and the energy gap /^E versus the 
impurity coupling strengths ai and a2 with double impurities at sites 1 and 4 for the two dimensional Ising 
lattice (7 = 1) in an external magnetic field h.—2. 

decays again as the impurity strength increases. The effect of ai in that case is less significant 
and makes the concurrence slowly decreases as ai increases which is expected since as the 
coupling between the two border sites 1 and 2 increases the entanglement between 1 and 4 
decreases. It is important to note that in general C(l,2) is much larger than C(l,4) since 
the border entanglement is always higher than the central one as the entanglement is shared 
by many sites. The entanglement between two regular sites is shown in fig. [TOT c) where the 
concurrence C(4,5) is depicted against ai and a2, the entanglement decays gradually as 02 
increases while ai has a very small effect on the entanglement, which slightly decreases as 
ai increases as shown. Interestingly, the behavior of the energy gap between the ground 
state and the first excited state of the Ising system AE versus the impurity strengths ai 
and ^2, which is explored in fig. ITW d) has a strong resemblance to that of the concurrence 



27 



0.2^ 




FIG. 21: (Color online) The concurrence C(l, 2), C(l, 4), C(4, 5) versus the impurity coupling strengths ai 
and a2 with double impurities at sites 1 and 2 for the two dimensional partially anisotropic lattice (7 — 0.5) 
in an external magnetic field h=2. 

C(4, 5) except that the decay of AE against 0:2 is more rapid. The effect of changing the 
location of the impurities is considered in fig. [SU] where the two impurities exist at sites 1 
and 4 in the Ising system. The behavior of the concurrences are very much the same as 
in the previous case except that the profiles of C(l,2) and C(l,4) have been exchanged as 
C(l,4) now represents the concurrence between the two impurity sites. 

The partially anisotropic system, 7 = 0.5, with double impurity at sites 1 and 2 and 
under the effect of the external magnetic field h = 2 is explored in fig. |2T1 As one can see, 
the overall behavior specially at the border values of the impurity strengths is the same 
as observed in the Ising case except that the concurrences suffer a local minimum within a 
small range of the impurity strength 0:2 between and 1 while corresponding to the whole ai 
range. The change of the entanglement around this local minimum takes a step-like profile 
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FIG. 22: (Color online) The concurrence C(l,2), C(l,4), C(4, 5) versus the impurity coupling strengths 
ai and a2 with double impurities at sites 1 and 2 for the two dimensional XY lattice (7 = 0) in an external 
magnetic field h=l. 

which is very clear in the case of the concurrence C(l,4) shown in fig. [2T](b). Remarkably, 
the local minima in the plotted concurrences coincide with the line of vanishing energy gap 
as shown in fig. [2T](d). which means that these minima correspond to a transition between 
a ground state and another one which takes place as the system parameters change. The 
anisotropic XY model with two impurities at sites 1 and 2 in an external magnetic field 
/i = 1.8 is explored in fig. [221 The entanglement for this system shows much sharper changes 
as a function of the impurity strengths and the sharp step changes take place in a narrow 
region of both ai and 02 specifically for — 1 < ai < 1 and —1 < 0:2 < 0. It is interesting to 
note that again the sharp step changes in the entanglement are corresponding to the line of 
minimum energy gap as shown in fig. [221(d) where this line varies continuously between a very 
small value and zero which explains the many steps appearing in the different concurrences 
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FIG. 23: (Color online) Dynamics of the concurrence C(l,2), C(l,4), C(5,7) with double impurities at 
sites 1 and 2 for the two dimensional Ising lattice (7 = 1) in an exponential magnetic field where a=l, b=2, 
and w=0.1. The straight lines represent the equilibrium concurrences corresponding to constant magnetic 
field /i = 2. The legend for all subfigures is as shown in subfigure (a). 

and particularly C(l,4) depicted in fig. [22r b). 
B. System dynamics with impurities 

Now we turn to the dynamics of the two dimensional spin system with double impurity 
under the effect of an external exponential magnetic field to test the ergodicity of the system 
as we vary the degree of anisotropy or the location of the impurities. In fig. [23] we consider 
the dynamics of the Ising system with two impurities at the sites 1 and 2 under the effect 
of an exponential magnetic field with parameters a = 1, 6 = 2 and a; = 0.1. As one 
can see the system shows an ergodic behavior for the different values of the impurities 
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FIG. 24: (Color online) Dynamics of the concurrence C(l,2), C(l,4), C(5,7) with double impurities at 
sites 1 and 2 for the two dimensional partially anisotropic lattice (7 = 0.5) in an exponential magnetic field 
where a=l, h—2, and w=0.f. The straight lines represent the equilibrium concurrences corresponding to 
constant magnetic field h = 2. The legend for all subfigures is as shown in subfigure (a). 

strengths ai, 0:2 = (0, 0), (0, 1), (1, 0) and (1, 1). The Ising system sustains its ergodicity for 
all shown values of impurities strengths and other tested values. The partially anisotropic 
system under the same condition behaves differently. Its pure case, ai, 02 = (0, 1) case and 
tti, 0.2 = (1, 1) case are all nonergodic as dipicted in figs. IMTa) andl2^b). Nevertheless, the 
system with impurity strengths ai = 1 and 02 = show ergodic behavior which means that 
the nonergodicity of the partially anisotropic system is sensitive for the strength and location 
of impurities. The isotropic system is explored in fig. [25] which behaves nonergodically for all 
impurity strengths. Testing the effect of the impurity location we consider the same system 
with impurities at the sites 1 and 4. While the Ising system shows ergodicty at all impurity 
strengths as shown in fig. [23l the partially and isotropic XY systems are nonergodic at the 
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FIG. 25: (Color online) Dynamics of the concurrence C(l,2), C(l,4), C(5,7) with double impurities at 
sites 1 and 2 for the two dimensional isotropic XY lattice (7 ~ 0) in an exponential magnetic field where 
a=l, b=:1.8, and w=0.1. The straight (thicker) lines represent the equilibrium concurrences corresponding 
to constant magnetic field h ~ 1.8. The legend for all subfigures is as shown in subfigure (a). 



different impurity strengths as plotted in figs. [27] and [28] respectively. 



V. CONCLUSION AND FUTURE DIRECTIONS 



We have investigated the nearest neighbor entanglement and ergodicity of a two- 
dimensional XY spin lattice in an external magnetic field h. The spins are coupled to 
each other through nearest neighbor exchange interaction J. The number of spins in the 
lattice are 7 where we may consider one or two of them as impurities. We have found that 
the completely anisotropic (the Ising), the partially anisotropic and isotropic systems behave 
in a very similar fashion to that of the one dimensional spin systems at the extreme, small 
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FIG. 26: (Color online) Dynamics of the concurrence C(l,2), C(l,4), C(5,7) with double impurities at 
sites 1 and 4 for the two dimensional Ising lattice (7 = 1) in an exponential magnetic field where a=f , b=2, 
and w^O.l. The straight lines represent the equilibrium concurrences corresponding to constant magnetic 
field h = 2. The legend for all subfigures is as shown in subfigure (a). 

and large, values of the parameter X = h/J but may deviate at the intermediate values. 
The first two systems show phase transition in the vicinity of the parameter critical value 
A = 2 and their entanglement vanishes as A increases. The entanglement of the isotropic 
system changes in a sharp step profile before suddenly vanishing in the vicinity of A = 2. 
The entanglement dynamics of the system with impurities was investigated under the effect 
of an external time-dependent magnetic field of exponential form. It was found that the er- 
godicity of the system can be tuned using the strength and location of the impurities as well 
as the degree of anisotropy of the coupling between the spins. It is interesting in future to 
investigate the same systems coupled to a dissipative environment and to examine the effect 
of impurity to tune the decoherence in the spin system and to investigate the ergodicity sta- 



33 



0.15 



CD 
O 

0) 

\ 

\ 

ID 
O 

o 
O 



0.1 



0.05 



(a) 



1'," ( ' ' 





10 



20 



30 



40 



50 



10 



20 



30 



40 



50 



-1> 



t (units of J"') 



t (units of J"^) 



0.2 

CD 

o 

S 0.15 



(c) 








a, = 1, a=l 




" 0.04 ^''"^^'>,v„ 
o 

o 

0.02 



'>v,Vw..i 



20 30 40 

t (units of J""") 



50 



0, 







10 



20 



30 



40 



50 



t (units Of J ') 

FIG. 27: (Color online) Dynamics of the concurrence C(l,2), C(l,4), C(5,7) with double impurities at 
sites 1 and 4 for the two dimensional partially anisotropic lattice (7 — 0.5) in an exponential magnetic field 
where a=l, h—2, and w=0.1. The straight lines represent the equilibrium concurrences corresponding to 
constant magnetic field h — 2. The legend for all subfigures is as shown in subfigure (a). 

tus under coupling to the environment. Furthermore we would like to investigate the same 
system with larger number of sites to test the system size effect and to clarify the critical 



value of the parameter A using finite size scaling |25l . |26| . Previously, the 19 sites triangular 
static Ising lattice was treated exactly using the the trace minimization algorithm 22] . The 
dynamics of entanglement in the 19-site XY system is currently under consideration, and 
by taking advantage of parallel computing we can reach 34 spins by far. 
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FIG. 28: (Color online) Dynamics of the concurrence C(l,2), C(l,4), C(5,7) with double impurities at 
sites 1 and 4 for the two dimensional isotropic XY lattice (7 = 0) in an exponential magnetic field where 
a=l, b=1.8, and w=0.1. The straight (thicker) lines represent the equilibrium concurrences corresponding 
to constant magnetic field h — 1.8. The legend for all subfigures is as shown in subfigure (a). 
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